Dynamical instability and dispersion management 
of an attractive condensate in an optical lattice 
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We investigate the stability of an attractive Bose-Einstein condensate in a moving ID optical 
lattice in the presence of transverse confinement. By means of a Bogoliubov linear stability analysis 
we find that the system is dynamically unstable for low quasimomenta and becomes stable near 
the band edge, in a specular fashion with respect to the repulsive case. For low interactions the 
instability occurs via long wavelength excitations that are not sufficient for spoiling the condensate 
coherence, producing instead an oscillating density pattern both in real and momentum space. This 
behaviour is illustrated by simulations for the expansion of the condensate in a moving lattice. 



The occurrence of dynamical instability and the possi- 
bility for dispersion management of Bose-Einstein con- 
densates (BECs) propagating in optical lattices have 
been the object of an intensive experimental and theo- 
retical activity in recent years [^. Dynamical instabil- 
ity, that originates from the interplay between nonlinear- 
ity of the interatomic interaction and periodicity of the 
lattice, has been proved to be a major mechanism for 
the disruption of the superfluid flow of repulsive BECs 
0, B i, II, i, 0, i, i, 111- Such a phenomena occurs 
via the exponential growth of arbitrarily small fluctua- 
tion when the relative velocity between the BEC and the 
lattice exceeds a critical threshold, above the critical ve- 
locity for Landau instability 3- the other side, in 
a regime of low interactions where nonlinear effects are 
negligible, it has been demonstrated that it is possible to 
control the dispersion of matter wave packets by tuning 
the velocity of the periodic lattice, in a similar fashion to 
the dispersion management in optics [ll|, [H, [l^ . This 
has striking effects for example on the expansion dynam- 
ics of the condensate, that can be reversed by switching 
from positive to negative effective mass regimes. 

The realization of BECs with tunable attractive inter- 
actions [fl [Tsl . Ill opens now interesting perspectives for 
the investigations of the effects of dynamical instability 
and dispersion management in unexplored regimes. 

In this paper we show that an attractive BEC moving 
in a ID optical lattice, contrarily to the repulsive case, 
is dynamically unstable in the first half of the Brillouin 
zone and can be stabilized near the band edge. Moreover, 
in case of weak nonlinearity and low quasimomenta the 
instability occurs via long wavelength excitations that 
do not destroy the condensate coherence, giving rise to 
density modulations that extend over several sites of the 
lattice. The effects of this interesting behaviour are dis- 
cussed by comparing the predictions of the Bogoliubov 
linear stability analysis with the solution of a time depen- 
dent Gross-Pitaevskii (CP) model [13] for the waveguide 
expansion of the condensate in the presence of the moving 
lattice. We show that in the unstable regime breathing- 
like oscillations are produced for low velocities, and that 
the absence of dynamical instability for higher velocities 



allows for a close inspection of negative mass effects near 
the band edge |12l] . 

Let us start by considering an infinite cylindrical con- 
densate of atoms of mass m, radially confined by the har- 
monic potential Vho{r) = (l/2)ma-'^r^, and subjected to 
the one-dimensional optical lattice V{z) = sEiiCOS^{kz). 
The parameter s represents the height of the lattice in 
units of the recoil energy = fi^k'^ /2m, k — 2-k jX be- 
ing the lattice wavenumber. The presence of the radial 
confinement makes convenient to use dimensionless ra- 
dial units. In the following, unless otherwise stated, we 
will express di stances in units of the harmonic oscillator 
length ar = \/'hlmLOr, time in units of lo~^ and energies 
in units of hiOr (keeping the same symbols in order to 
simplify the notation). 

The behaviour of the system is governed by the 3D 
GP equation for the condensate wavefunction [ITj] . How- 
ever, owing to the radial confinement it can be conve- 
niently described also by means of an effective ID model, 
the Non-Polynomial Schrodinger Equation (NPSE) [T^ . 
Though onedimensional, this model includes an effective 
radial-to-axial coupling, and has been shown to pro vide 
a realistic description in several situations [1, H, Il8l | . It 
is obtained from the GP equation by means of a factor- 
ization of condensate wavefunction in the product of a 
Gaussian radial component of z- and t-dependent width, 
(j{z, t), and of an axial wave function t) solution of 
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with (T^(z,t) = \/l + g\f{z, The coupling constant 
g is defined as g = 2Na/ar, a being the scattering length 
and N the number of atoms in each site of the lattice. 
The wavefunction (p is normalized to unity over a single 
period. In the limit of vanishing lattice, s = 0, iy9 is spa- 
tially uniform and the normalization yields \(p\'^ — 2/A. 

Equation ([1]) admits stationary solutions in the form 
of Bloch waves (p{z,t) = 0p(z)e**^''*~^^) where (j)p has the 
same period of the lattice, the Bloch wave vector p rep- 
resents the quasimomentum of the condensate and is re- 
stricted to the first Brillouin zone, and ^ is the chemical 
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FIG. 1: (Color online) Critical values \gcr\ of the interaction 
strength as a function of s, for three values of A. For each 
value of A the system is stable against collapse in the region 
below the corresponding line for \gcr\- 



potential. The stationary NPSE for (^^ becomes 



(2) 



with i/p = (idz — p)^/2 + si?fl cos^(7rz/A). The behavior 
of the system is therefore characterized by three dimen- 
sionless parameters: A, s and g. 

In case of attractive interaction among the atoms, 
g — — |f/|, the interaction strength is limited by a crit- 
ical value gcr at which the system collapses. The value of 
gcr depends on both A and s, and its behaviour is shown 
in Fig. [TJ For s = the stability is determined by the 
presence of the radial trapping, and the critical value is 
ffcrl — y^/'i.. For higher values of \gcr\ the energy func- 
tional becomes unbounded from below and the equation 
for no longer admits a real solution. By increasing 
s also the effect of the lattice becomes important and 
reduces the critical value of \g\. 

Though stable against collapse the above stationary 
solution can present other kind of instabilities, among 
which energetic (or Landau) and dynamical instability 
(DI). The latter can be investigated by means of a lin- 
ear analysis, considering time dependent perturbations of 
the form ip{z,t) = e'C/^*-?^) [0p(^) %(z, i)] where the 
fluctuations 5(j)p can be expanded in terms of Bogoliubov 
quasiparticle amplitudes Upq,n{z) and Vpq,n{z) with the 
same periodicity of the lattice 
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q being the quasimomentum of the excitations and n the 
Bloch band index. The Upq^n and Vpq^n are solutions of 
the Bogoliubov-like equations that corresponds the 
eigenvalue problem for the 2x2 matrix M defined as 

M = {A+,B; -B*, -A_), with [9] 
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DI occurs when the frequencies lOpq of some modes get 
a nonzero imaginary part, indicating an exponential 
growth. If initially occupied, these modes rapidly become 
macroscopically populated driving the condensate away 
from the unperturbed state [1, d, H, Q . In a similar way, 
but considering the contribution bE to the energy func- 
tional of time independent perturbations, one can define 
energetic instability, occurring when a^^M has negative 
eigenvalues (0-3 being the third Pauli matrix) [1, Q . 

The Bogoliubov frequencies iOpq can be obtained by 
numerical diagonalization of the matrix M and the cor- 
responding dynamically unstable regions can be drawn 
in the p-g plane where the imaginary component of LOpq 
is nonvanishing. Since complex frequencies always come 
in conjugate pairs, the appearance of an imaginary com- 
ponent corresponds to a resonance between the phonon 
and antiphonon branches of the spectrum 0. 

Let us now focus to the case of an attractive conden- 
sate. To illustrate a specific case we take A = 2.5, but 
the efl:ects we are going to discuss can be obtained in a 
wide range of parameters, with different values of the ra- 
dial confinement and lattice periodicity. Typical results 
for the stability diagrams are shown in Fig. [5] for differ- 
ent values of s and g. In this picture we also show the 
regimes of energetic instability obtained from the diag- 
onalization of (T3M (shaded regions), indicating that in 
the attractive case the system is energetically unstable 
for any velocity, and therefore it is not a superfluid ac- 
cording to the Landau criterion (l7| . 

As regards DI, these diagrams show a remarkable be- 
haviour compared to the case of repulsive interactions 
investigated so far 0, Q (see e.g. Fig. [2^-i), and there 
are at least three points that make the attractive case 
particularly interesting, i) Increasing the lattice velocity 
stabilizes the system: attractive condensates are dynami- 
cally unstable at low quasimomenta, but can be stabilized 
above a certain critical velocity. This is just the oppo- 
site behaviour of that for repulsive BECs, that are stable 
below a critical velocity 2, 8]. ii) For low quasimomenta 
and shallow lattices DI takes place via long wavelength 
(low q) excitations, see Fig. [5J In this regime the excita- 
tions spectrum has finite real and imaginary components 
and reads LOpq ~ qa(j)) + i(3{q) (the origin of this form will 
be clear from the discussion of the s — case here below) . 
These modes are not sufficient to produce a site-to-site 
dephasing, but are rather expected to induce collective 
oscillation of the system, as will be shown later, iii) Neg- 
ative effective mass and DI appear in separate regions 
of the p-q plane, thus opening the possibility for being 
investigated in a more clear way. For example, one can 
tune the dispersion relation to negative values by loading 
the condensate in a quasimomentum state above (73/2, 
avoiding the disturbance of DI that prevent a close in- 
spection near the band edge in the repulsive case d, . 

Following we can compare these numerical results 
with the analytics in two special cases: vanishing lattice 
(s = 0) and noninteracting limit {g = 0). For s = the 
diagonalization of the matrix M yields the following an- 
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FIG. 2: (Color online) Stability diagrams as a function of the 
quasimomenta p and q of the condensate and the excitations, 
for different values of s and g. Colored regions indicate where 
the system is dynamically unstable. The colour scale corre- 
sponds to the growth rate of the unstable modes for a given 
p, q (imaginary part of cupq). The solid and dashed lines in (a), 
(b) and (c) correspond to the analytic predictions for s = 0, 
whereas the bold line in (g), (h) and (i) refer to the case g — 
(see text). Light shaded area: energetic instability. Dotted 
region in (g), (h), and (i): regimes of DI for the repulsive case 
9 = \9\- 



alytic expression for the excitation spectrum, lu± ^ pq± 
^/^^JIT^, with c = (3(tV2 - 1/0-2 - i/2cr2)/2. This 
expression for lu± is formally the same of the strictly ID 
case , but with a functional dependence of the interac- 
tion term c on the condensate density (through a). 
In case of attractive interactions, g < 0, it is straightfor- 
ward to verify that c is always negative. In this case com- 
plex frequencies appear for < 2\c\ (for any p) and for 

P = W{q-2qBY/\ + c{q-2qBY - ^ q^ / A + cq^) /2qB, 
as indicated in Fig. [2K-c by dashed and solid lines re- 
spectively. In the former case, that is peculiar of attrac- 
tive interactions, both real and imaginary component of 
bjpq are finite, whereas in the latter the frequencies are 
purely imaginary (as for repulsive condensates 0, Q)- It 
is easy to verify that the former condition (g^ < 2|c|) 
corresponds also to the the presence of negative eigen- 
value of (73 M even at p = 0, the signature of a Landau 
instability. 

In the noninteracting limit c = the excitation en- 
ergies can be directly obtained from the noninteracting 
spectrum of the condensate E{p), u!±{q) = ±{E{p ± 
q) — E{p)). In this case the resonance condition reads 
E(j) + q) - E{p) = E{p) - E{p - q) [1]. The points 
that satisfy the above constraint are shown by the line 
in Fig. and nicely reproduce the boundary be- 

tween the unstable regions for the attractive and repul- 
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FIG. 3: (Color online) Density plot of the condensate wave- 
function along the axial direction (top) and of its momentum 
distribution (bottom) as a function of time during the expan- 
sion in the waveguide, in the presence of a lattice moving at 
V = 0.2vb (see text). 



sive cases. For low q's this resonance condition can be 
related to the condition for the vanishing of the inverse 
effective mass 1/m*: in this case it can be easily demon- 
strated that besides the trivial solution Ep — Ep, one also 
has g2 ^ 12£;(2)(p)/^(4)(p) (where £;(")(p) = d''E/dp'') 
leading to 1/m* = in the limit q — 0. This result pro- 
vides an explanation of the numerical coincidence found 
in concerning this point. 

To illustrate how the features of the spectrum dis- 
cussed so far affect the system dynamics we now discuss 
the expansion of the condensate in a waveguide in the 
presence of a moving lattice. As an example we con- 
sider an elongated condensate composed by iV = 1000 
atoms of '^^K with attractive interactions characterized 
by a scattering length a = — ao, tuned by means of Fes- 
hbach resonances [l6|. The BEC is initially prepared in 
a waveguide of radial frequency LUr = 27r x 1.5 kHz, in 
the presence of an additional harmonic confinement of 
frequency lu^ = 2tt x 20 Hz along the waveguide. 

Then, at i = 0, the axial confinement is released and 
the condensate is let expand in the waveguide while a op- 
tical lattice of wavelength A = 1040 nm, moving at veloc- 
ity V, is ramped up to a final intensity s = 5 that is kept 
constant for the subsequent dynamics. As ramp time we 



choose tramp = 10 ms, that is sufficient to project adi- 
abatically the condensate in a state of quasimomentum 
p/pb v/vb [S, 12]. This scheme is similar to experi- 
mental setup in [l^, with the difference that there also 
the radial confinement was removed in order to reach 
rapidly the non-interacting regime. Conversely, since 
here we want to probe the effects of nonlinearity, the 
presence of the waveguide helps to keep interactions for 
longer times. 

The dynamics is investigated by solving the time de- 
pendent NPSE in Eq. ^ for different velocities of the 
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FIG. 4: (Color online) (a) Evolution of the axial rms size 
of the condensate in the regime of negative effective mass, 
for different velocities, compared to the free expansion in the 
waveguide (solid line). The condensate is let expand freely 
for 10 ms in the waveguide before ramping on the lattice, (b) 
Effective mass vs quasimomentum. (c) Axial density plot of 
the expanding BEC for v = O.Svb- 



lattice [l9[. Let us first consider the case v = 0.2vb 
shown in Fig. [3l where we plot the evolution of the ax- 
ial density and the momentum distribution in the refer- 
ence frame co-moving with the condensate at the end of 
the ramp. According to the above linear stability analy- 
sis we expect the system to be dynamically unstable via 
long wavelength phonons. In particular, since the con- 
densate occupies about 30 sites with an average number 
of atoms per site N « 30, the stability behaviour should 
be close to that shown in Fig. [2^. Indeed the simula- 
tions show a rapid population of modes at small quasi- 
momenta around q = O.OSSgs, corresponding to density 
modulations over several sites of the lattice, and match- 
ing the most unstable modes of the uniform system in 
Fig. [2^. Remarkably, the growth of these modes does 
not induce decoherence as observed so far with repul- 
sive BECs j8|l, but produces instead an oscillatory be- 
haviour that is maintained for several periods. These 
oscillation can be accounted for by the finite real com- 
ponent of the excitations spectrum, and by taking into 
account the momentum spread Ap due to finite size of the 
system. Indeed, the fitted frequency v = 44.5 Hz in Fig. 
13] is consistent with the real part of frequency of the most 
unstable modes at q. Such a behaviour is conserved until 
the lattice velocity is tuned close to w ~ 0.5vb, where the 
presence of unstable modes with q ~ qs favor a change 
of the condensate phase between neighboring sites of the 
lattice, leading to a complex dynamics that rapidly spoils 
the phase coherence of the BEC. 

Another interesting scenario shows up in the dynami- 
cally stable region for q > 0.5qB where the the effective 
mass m* changes sign. As discussed and experimentally 



demonstrated in [1^, |l^ a change of sign of m* corre- 
sponds to a time-reversed evolution and is therefore ex- 
pected to produce a contraction of the axial size of a 
condensate initially expanding outwards TT]. This is in- 
deed shown in Fig. [H^a) where we plot the evolution of 
the axial rms size for different velocities. This picture 
shows that in the presence of the lattice the expansion 
is slowed down until the wave packet eventually reaches 
its minimum allowed size Zmim and then the dynamics 
inverts again. This effect is more evident for lower veloc- 
ities where a larger negative m* (see Fig. IHb)) produces 
a hampering of the expansion for a longer time. Con- 
versely, the expansion is enhanced near the band edge 
due to the lower value of \m*\. Notably this enhance- 
ment is favored also by the attractive interactions that 
are turned into an effective repulsion when m* changes 
sign, causing a boost of the expansion with respect to 
the free case (even when m < \m*\). This is the same 
mechanism that limits the value of Zmm compared the 
noninteracting case [ij]. Finally, Fig. dfc) demonstrates 
that in this regime of negative effective mass the system 
is not affected by DI, as discussed above. 

In conclusion, we have discussed the occurrence of dy- 
namical instability for attractive BECs moving in a ID 
optical lattice in the presence of transverse confinement 
finding a remarkable behaviour with respect the usual 
repulsive case. Contrarily to the latter, attractive con- 
densates are dynamically unstable at low quasimomenta 
and can be stabilized above a certain critical velocity. 
By suitably adjusting the lattice height and velocity it is 
also possible to tune the disrupting character of dynam- 
ical instability [1] into a breathing-like collective oscilla- 
tion of the system induced by long wavelength phonons. 
Moreover, since negative effective mass and dynamical 
instability appear in separate regions of the p-q plane of 
the condensate and excitations quasimomenta, this al- 
lows for a more clear investigation of both effects. We 
have also explicitly verified that the results presented 
here by means of the NPSE model hold even in a full 3D 
CPE approach, provided that the radial confinement is 
tight enough to prevent the mixing of the radial branches 
of the spectrum (as the case considered here) [§]. These 
results may be particularly relevant for current experi- 
ments with Bose-Einstein condensate with interactions 
tunable around zero 111, UK M • 
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